Vortex motion in a finite— size easy— plane ferromagnet and application to a nanodot 
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We study the motion of a non-planar vortex in a circular easy-plane ferromagnet, which imitates 
a magnetic nanodot. Analysis was done using numerical simulations and a new collective variable 
theory which includes the coupling of Goldstone-like mode with the vortex center. Without magnetic 
field the vortex follows a spiral orbit which we calculate. When a rotating in-plane magnetic field is 
included, the vortex tends to a stable limit cycle which exists in a significant range of field amplitude 
B and frequency ui for a given system size L. For a fixed w, the radius R of the orbital motion is 
proportional to L while the orbital frequency £1 varies as 1/L and is significantly smaller than u. 
Since the limit cycle is caused by the interplay between the magnetization and the vortex motion, 
the internal mode is essential in the collective variable theory which then gives the correct estimate 
and dependency for the orbit radius R ~ BL/uj. Using this simple theory we indicate how an ac 
magnetic field can be used to control vortices observed in real magnetic nanodots. 

PACS numbers: 75.10.Hk, 75.30.Ds, 05.45.-a 
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I. INTRODUCTION 

Nonlinear topological excitations in 2D spin systems 
of soliton or vortex type are known to play an essential 
role in 2D magnetism. For example, solitons break the 
long-range order in 2D isotropic magnets. Vortices play 
a similar role in 2D easy-plane magnets. Magnetic vor- 
tices have been studied since the 1980s. They are impor- 
tant for the dynamical and thcrmodynamical properties 
of magnets, for a review see Ref. OJ The vortex contri- 
bution to the response functions of 2D magnets has been 
predicted theoretically 2 and observed experimentally 3 . 

A second wind in the physics of magnetic vortices ap- 
peared less than five years ago due to the direct observa- 
tion of vortices in permalloy (Py, Ni 8 oFe2o) 4 ' 5 i 6 i 7 i 8 i 9 and 
C o 10 ' 11 ! 12 magnetic nanodots. Such nanodots are sub- 
micron disk-shaped particles, which have a single vor- 
tex in the ground state due to the competition between 
exchange and magnetic dipole-dipole interaction^ A 
vortex state is obtained in nanodots that are larger 
than a single domain whose size is a few nanometers 
(e.g. for the Py-nanodot the exchange length Z cx = 
5.9 nm). The vortex state of magnetic nanodots has 
drawn much attention because it could be used for high- 
density magnetic storage and miniature sensors^i^. For 
this one needs to control magnetization reversal, a pro- 
cess where vortices play a big rofoi&. The vortex signa- 
ture has been probed by Lorentz transmission electron 
microscopyii*ii and magnetic force measurementsiS*^. 
Great progress has been made recently with the possi- 
bility to observe high frequency dynamical properties of 
the vortex state magnetic dots by Brillouin light scatter- 
ing of spin wavesi 9 ^, time-resolved Kerr microscopy 9 , 
phase sensitive Fourier transformation technique^!, and 



X-ray imaging technique^. These have shown that the 
vortex performs a gyrotropic precession when it is ini- 
tially displaced from the center of the dot, e.g. by an 
in-plane magnetic field pulseASiSi 

In general the vortex mesoscopic dynamics is described 
by the Thiele collective coordinate approacb J 2 ^Si 2 i22i, 
which considers the vortex as a rigid structure not hav- 
ing internal degrees of freedom. 1 However recent experi- 
mental and theoretical gtudiesiiiaS^iSSiSLi&iSSi indicate 
phenomena which can not be explained using such a sim- 
ple picture. One striking example is the switching of the 
vortex polarizatiorAiii 3 ^ 3 ^ 3 ^ 3 ^, where coupling occurs 
between the vortex motion and oscillations of its core. 
Another one is the cycloidal oscillations of the vortex 
around its mean pathS 9 ^^ where the dynamics of the vor- 
tex center is strongly coupled to spin waves. In this way 
the internal dynamics of the vortex plays a vital part. 
One of the first attempts to take into account the inter- 
nal structure of vortices was presented in Ref. |35| which 
showed that a variation of the core radius slaved to the 
position explained the motion of a vortex pair across an 
interface between two materials of different anisotropy. 
Some progress has been achieved in Ref. |3fj where we 
have confirmed that internal degrees of freedom play a 
crucial role in the dynamics of vortices driven by an ex- 
ternal time-dependent magnetic field in a classical spin 
system. 

Here we present a complete study of this problem using 
direct numerical simulations of the spin system and a col- 
lective variable theory which includes an internal mode. 
We show that the periodic forcing of the system by the 
time dependent magnetic field together with the damp- 
ing stabilizes the vortex in a finite domain. This limit 
cycle exists because of the interplay between the magne- 
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tization and the vortex position so that it is essential to 
include an internal mode in the collective variable theory 
to describe it. When this is done, the theory yields the 
domain of stability in parameter space and the main de- 
pendencies on the field amplitude B and frequency to. It 
can be seen as a one of the first generalizations to vortices 
of the collective variable theories developed for ID Klein- 
Gordon kinks by Rice^LSS^S which include the width of 
the kink together with its position. 

In the next section [n] we formulate the continuum 
model, discuss the role of different types of interactions 
and briefly review the main results on the structure of 
the vortex solution. The vortex motion without external 
field is examined in section IIIII It follows a spiral or- 
bit as a result of the competition between the gyroforce, 
the Coulomb force and the damping force. In section 
^3 with the ac driving, numerical simulations show that 
the vortex converges to a stable limit cycle. We give its 
boundaries in parameter space and indicate how the ra- 
dius and frequency of the vortex orbital motion depends 
on the field and geometry parameters. SectionlVlpresents 
and discusses in detail the new collective variable theory 
of the observed vortex dynamics which takes into account 
the coupling between an internal shape mode and the 
translational motion of the vortex position. In section 
IVII we link this with the individual spin motion observed 
in the simulations and indicate how these effects can be 
observed and used in real nano magnets. 

The model we consider is a ferromagnetic system with 
spatially homogeneous uniaxial anisotropy, described by 
the classical Heisenberg Hamiltonian 



J 
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(S Z r , 



(1) 



Here S n = (5^,5^,5*) is a classical spin vector with 
fixed length S on the site n of a two-dimensional square 
lattice, and the exchange integral J > for a ferromag- 
net. The first summation runs over nearest-neighbor 
pairs (n, n'). We assume a small anisotropy leading to an 
easy-plane ground state. This anisotropy can be either 
of the exchange type, with < 6 <C 1, or of the on-site 
type, with < K < J. 

Extending ideas of Ref . |3^ we study the movement of a 
vortex in this system under the action of a magnetic field 
B{t) — (B cos tot, -Bsinw£,0), which is spatially homoge- 
neous and is rotating in the plane of the lattice. This 
field adds an interaction of the form 



f(t) = -jB 0S£ cos "t + ft sin tot) 



(2) 



where 7 = 2/is/fi is the gyromagnetic ratio. 

The spin dynamics is described by the Landau-Lifshitz 
equations with Gilbert damping 



dS n 
~dt 



dS r , 



dt 



(3) 



where J4? = Jf +-f(t) is the total Hamiltonian. Eqs. Q 
preserve the length of the spins | S n \ = S, which has units 
of action. Another form of Eqs. © more suitable for spin 
dynamics simulations is given in Appendix lAl 



II. CONTINUUM LIMIT 

In the case of weak anisotropics 5 <C 1, K <C J, t he 
characteristic size of excitations Iq — ay/ J/{AJ8 + K) is 
larger than the lattice constant a, so that in the low- 
est approximation on the small parameter a/lg and weak 
gradients of magnetization we can use the continuum ap- 
proximation for the Hamiltonian |JT|I 



TC2 

H = J%-Eo = — I d 2 x 



(V*) 



I 2 
'0 



(4) 



where Eq is a constant. The spin length has been rcscalcd 
so that 



5/5= (VI 



^;VT 



(5) 



is a unit vector. The length Iq coincides with the radius of 
the vortex core obtained in Ref. |2^ for on-site anisotropy 
type alone (5 = 0). For the case of exchange anisotropy 
alone (K = 0), it is also customary to use the length 
r v = ay/ (1 — 6)/46, which is obtained from an asymp- 
totic analysis and is to be identified later with the radius 
of the "core" of a vortex4fi^i However, for the range of 
5 we are interested in, i.e. for 6 < 0.1, the difference 
between r v and la is negligible. 

The interaction with a homogeneous time-dependent 
magnetic field is expressed as 



V(t) 



JS 2 

I 2 
'() 



d 2 x{b(t) ■ s(r,t)) 



= -JS 2 b J d 2 £ Vl-m 2 cos - vt) 



(G) 



In order to simplify notations we use here and be- 
low the dimensionless coordinate £ = r/lo, the dimen- 
sionless time r = toot, the dimensionless driving fre- 
quency v = uj/luq and the dimensionless magnetic field 
6 = 7 B/w r2i^where 



uj a = S (AJ5 + K) 



(T) 



In all real magnets there is, in addition to short-ranged 
interactions, a long-ranged dipole-dipole interaction. In 
the continuum limit this interaction can be taken into 
account as energy of an effective demagnetization field, 



d 2 xM-H (m \ 



where M is the magnetization. Generally, this field is 
a complicated functional of M . However, in the case of 
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a thin magnetic film (or particle) the volume contribu- 
tion to the demagnetization field is negligible, and only 
surface fields are important. The face surfaces produce 
a local field H^ m ' = — 47rA/o e 2 for the sample with the 
saturation magnetization Mq. Then the dipole-dipole 
interaction can be taken into account by a simple re- 
definition of the anisotropy constants, K — > K cS = 
K + 4irM§a 2 / S 2 , leading to a new magnetic lengtbii 



Without magnetic field the ground state of the system 
is a uniform planar state m = and 4> — const. The 
field changes essentially the picture: spins start to precess 
homogeneously in the XY-plane, cj> = ip + vt. Such a 
precession causes the appearance of a z-component of 
magnetization, m — const. From Eqs. (|13|) . we find that 
the equilibrium values of m and <f> satisfy the following 
equations, 



/ ; cff - 



J 



4J6 + K + 4iTM$a 2 /S 2 ' 



(8) 



This is the case of so-called configurational or shape 



anisotropy 



5.14.45 



The lateral surface affects only the 



boundary conditions, see Refs. Hi 17 for details. For ex- 
ample, for a very thin magnetic particle, which corre- 
sponds to our 2D system, free boundary conditions are 
valid, and we will use them in the paper. 

Thus, the total energy functional, normalized by JS 2 , 
reads 



I 



d 2 Z 



m 
T 



(b-s) 



(9) 



where we have rescaled the magnetic length in accordance 
with |(|J. 

The continuum version of the Landau-Lifshitz Eqs. 
becomes 



dcf> 5$ e dm 

dr 5m (1 — m 2 ) dr ' 
dm 8$ 2 d<j> 



(10a) 
(10b) 



v 
m 
b sine 



e 2 v 2 



ev 



1 9 " 

1 — m z 
^2 = 0, 



(14a) 
(14b) 



so that this state can only exist if b > ev (otherwise 
only the ground state with m = and cp — const exists) . 
Assuming m <C 1, we obtain 



1 - VP - 



£ 2 V 2 



vt + 7r + arcsin ■ 



ev 



VT 



(15) 

Note how the magnetization m is proportional to the 
field frequency v so that its sign is important. Below 
we discuss the role of this homogeneous solution in the 
vortex dynamics. 

The continuum analogue of the power-dissipation rela- 
tion (|A2|) for the total energy functional <S[s] is calculated 
from Eqs. lfTT|l . ltT2*|l and gives 



dS 
~dV 



db 



(16) 



Formally Eqs. (|14f) have two solutions. One can check 
that only for the solution (|15fl the dissipation balances 
the work done by the field, so that the energy S tends to 
be stabilized. 



These equations can be derived from the Lagrangian 

/ : : / d 2 i{\-m) d -^--S\ S \ (11) 



and the dissipation function 



d 2 ^ 



(12) 



1 — m 2 



dm 
~dr~ 



+ (1 - m 2 ) 



Then Eqs. I|10(l result explicitly in 



m (Vm) 2 
^ — + rr 

(1 — m 2 ) 2 

bm cos(0 — vt) 
\/l — m 2 



7h 



Am 



(1 — m 2 ) 8t 



1 — m 2 

dm ho ^ 

— » ( 13a ) 



dm 



— = V[(l - m 2 ) V<f>] - bVl-m 2 sin(0 - vt) 



e (l_ m ^)_ 



(13b) 



Static Vortices 

The simplest nonlinear excitation of the system is 
the well-known non-planar magnetic vortex. We re- 
call briefly the structure of a single static vortex at zero 
field. In this case the pair of functions (m, 4>) satisfies the 
Eqs. I|10|) with the time derivatives set to zero and 6 = 0. 
If we look for planar solutions (m = 0) for the <fi field, 
Eq. I|10bl) becomes the Laplace equation. For the vortex 
solution located at Z = X + iY = i?exp(i$) the 0-field 
has the form: 



= ipo + q arg (z - Z) . 



(17) 



where z = x + iy is a point of the XY-plane, q G Z is 
the 7Ti topological charge of the vortex (vorticity). We 
will call the solution with q = 1 a vortex and the solution 
with q = — 1 an antivortex. The expression (|17JI does not 
satisfy the boundary conditions for a finite system. For 
our circular system of radius L (in units of Iq) and free 
boundary conditions the solution is^4 

<j) = arg (z - Z) - arg (z - Z T ) + arg Z, (18) 
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where the "image" vortex is added at Zi = ZL 2 /R 2 to 
satisfy the Neuman boundary conditions. The last term 
in 1)18(1 is inserted to have the correct limit for L — > oo. 

The m-field has radial symmetry, m = 
cosft (p = \z — Z\). From l(13a() and (fTTjl one can 
derive that #(•) satisfies the following differential 
problem: 



d 2 e 1 do 

dp 2 p dp 
cos 6(0) = p. 



sin ft cos ( 



1 - 



cos#(oo) 



1 

7 

0. 



= 0, 



(19a) 
(19b) 



where p = ±1 is the so-called polarity of the vortex. 
The solution of this differential problem is a bell-shaped 
structure with a width in the order of ^o- 



force appears as a result of boundary conditions, it de- 
scribes the 2D Coulomb interaction between the vortex 
and its image 



S 



inl 



(oq + 7T In 



L 2 -R 2 



(23) 



where S'q ~ ir is the energy of the vortex core^ 

In order to generalize the effective equations of the 
vortex motion for the case of the magnetic field we de- 
rive now the same effective equations by the effective La- 
grangian technique as it was proposed in Refs. |.".:|:; r i:il 
Inserting Ansatz 1)20(1 into the "microscopic" Lagrangian 
((11(1 and the dissipative function ((12(1 . and calculating 
the integrals, we derive an effective Lagrangian (see Ap- 
pendix El for the details) 
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III. VORTEX MOTION AT ZERO FIELD 

A standard description for the steady movement 
of magnetic excitations was given first by ThieleiSSi2& 
Huber— , Nikiforov and Sonin^ have first applied this 
approach to the dynamics of magnetic vortices, using a 
traveling wave Ansatz s(z,t) — s (z — Z(t)). In terms 
of the fields m and 4> such an Ansatz is 

m(z,r) = cosft (|z- Z(r)\) , (20a) 
4>(z,t) = arg(z - Z{t)) - arg (z - Z/(r)) + argZ(r), 

(20b) 

where the function 6*(») describes the out-of-plane struc- 
ture of the static vortex, and is the solution of Eqs. ((19(1 . 

To derive an effective equation of the vortex motion 
for the collective variable R(t) = (X(t), Y(t)), we 
project the Landau-Lifshitz Eqs. 1(10(1 over the lattice 
using Ansatz ((20(1 . We obtain a Thiele equation in the 



G e 2 x R - 2itt]R + F = 0, 



(21) 



where the dot indicates derivative with respect to the 
rescaled time r. The first term, the gyroscopic force, acts 
on the moving vortex and determines the main properties 
of the vortex dynamics. The value of the gyroconstant is 
well-known G = 2irpq^^ in our case for the vortex with 
positive polarity and unit vorticity G — 2-k. The second 
term describes the damping force with a coefficient 27 ' 48 



r) = ie(lnL + Ci), 



(22) 



where C\ ~ 2.31 is a constant coming from the m field 
and is calculated in the appendix, see formula ((B13(l . The 
In L dependence in 7/ was obtained in Ref . 1271 

The last term in Eq. 1(21(1 is an external force, acting 
on the vortex, F = —^rS', where $ is the total energy 
functional @. Without magnetic field (b = 0) such a 



(24) 



In the same way we derive the effective dissipative func- 
tion 



& = irr/R 2 = ttj] (h 2 + R 2 <t> 2 



(25) 



The equations of motion are then obtained from the 
Euler-Lagrange equations 



dJ? _d_ (d^\ _ 



for the Xi = {R,&}, 



1 



• R_ 

+ V R~ L 2 - i? 2 ' 

R A 
R - 7 ^ 



(26) 

(27a) 
(27b) 



This set of equations is equivalent to the Thiele Eq. 121(1 , 
when going to polar coordinates. 

For zero damping (e = 0) two radial forces act on the 
vortex (gyroforce and Coulomb force) and compensate 
each other, providing pure circular motion of the vortex. 
In that case the radius R of the orbit is arbitrary. Using 
Eqs. 1(27(1 it is easy to calculate the frequency of this 
circular motion for a given R, see Ref. 01 



Q{R) = 



1 



L 2 ~R 2 ' 



(28) 



When the damping is present, there appears an addi- 
tional damping force which can not be compensated by 
other forces. Thus the trajectory of the vortex becomes 
open-ended, following the logarithmic spiral from l(27b(l : 



$-$ = - In— , 
V R a 

where Rq and $o are constants. 



(29) 
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FIG. 1: Trajectory of a vortex at zero field for a time interval 
< t < 10 4 . The damping e = 0.01 was switched off at 
t — 1600. The vortex with q = p = 1 was launched from 
Z — 20.5a + i23.5a on a lattice of radius L = 20a ps 11.3. 
"Clean" circular trajectories, where the vortex is free of spin 
waves, are obtained with this method. In the whole study the 
anisotropy is set to 8 = 0.08. The damping is e = 0.01. 



IV. NUMERICAL SIMULATIONS OF THE 
VORTEX DYNAMICS 

To investigate the vortex dynamics, we integrate nu- 
merically the discrete Landau-Lifshitz equations (|A1|) 
over square lattices of size (2L) 2 using a 4th-order 
Runge-Kutta scheme with time step 0.01. Each lattice 
is bounded by a circle of radius L on which the spins are 
free corresponding to a Neuman boundary condition in 
the continuum limit. In all cases the vortex is started 
near the center of the domain and the field and damping 
are turned on adiabatically over a time interval of about 
100. We have only considered vortices of fixed polarity 
p = 1. More details on the numerical procedure and in 
particular the vortex tracking algorithm can be found in 
Ref.Hl 

We have fixed the exchange constant J = 1 as well 
as the spin length 5=1. All cases presented here are 
for the anisotropy 5 — 0.08, corresponding to Iq rs 1.77a 
so that we are close to the continuum limit. The lattice 
radii we consider here are 20a < L < 100a. 

To validate the simple theory presented in the previous 
section we considered the case with no magnetic field. In 
the absence of damping the vortex should follow a circu- 
lar orbit and its frequency of rotation should be given by 
(|28|1 . Starting with a vortex initial condition for m and <fi 
given by l|20|l . it is possible to "prepare" circular trajec- 
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FIG. 2: (Color online) Azimuthal angle 3> of the vortex po- 
sition as a function of the logarithm of the radial position R 
for four different values of the damping e. The lattice radius 
is L = 20a w 11.3. 



tories of arbitrary radius by applying damping. This kills 
all spin waves coming from the imperfect initial condition 
and drives the vortex to the selected radius following the 
spiral (|29|l . Once the chosen radius is reached, damping 
is turned off adiabatically over a time greater than 100 
(1/e) and the vortex will keep its circular orbit indefi- 
nitely Such a scenario is shown in Fig. 

We now analyze the spiral trajectories obtained when 
damping is present. In Fig.|21we plot the measured angle 
of rotation $ (in radians) as a function the logarithm 
of the measured radius R for four values of damping. 
The vortex is started every time from the same place 
($0 = 7r/2, Ro = 3a ) in the lattice. The behavior given 
by the spin simulation shown by full lines agrees well 
with the relation l|29|) given by dashed lines. Note that 
the constant C\ is important to obtain a quantitative 
agreement because it is of the same order as the term 
InL. 

To study the vortex dynamics in the presence of the 
rota ting field, we extend the simulations described in 
Ref. y3|. There we investigated the dynamics of the out- 
of-plane structure of the vortex, focusing on the phe- 
nomenon of switching, which occurs when vp < 0. Here 
we consider vortices with positive polarity p — 1 and 
v > so that no switching occurs. 

For simplicity we fixed the damping e = 0.01 in (J3J) 
and varied the parameters (b, v : L). We checked that the 
effects reported here occur for a range of anisotropics 
and damping around these values. Given a combination 
of the parameters {v,b) of the field, the radius L of the 
system and the damping e, we have observed that cither 
the vortex escapes from the system through the border 
or it stays inside for all times. In the latter case, it can 
approach a limit cycle for a broad range of the field pa- 
rameters. Fig. 13 shows two vortex trajectories starting 
from different positions and converging to the same cir- 
cle. When the limit cycle exists, its basin of attraction is 
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FIG. 3: (Color online) Two trajectories of a vortex from 
simulations of the many-spin model GJ-Q, on a lattice of 
radius L — 78a » 44, with a rotating field (y — 0.125, b — 
0.002). For this field, all trajectories converge to the same 
circle independently of the vortex's initial position, provided 
it is not too close to the system border. 
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FIG. 4: (Color online) Diagram of types of trajectories in 
the (to, B) parameter plane corresponding in the (y, b) plane 
to the range < v < 0.3 and < b < 0.033. The radius of 
the system is L — 36a m 20. The term "confined" means that 
no limit cycle was reached, though the vortex stayed inside 
the system, during the time of observation r < 6400. 



very large as can be seen by starting the vortex at differ- 
ent positions and seeing it converge to the same circle. In 
other words, the system keeps no memory of the initial 
position of the vortex. 

To exist, the limit cycle needs both magnetic field and 
damping: once it is attained, switching off or changing 



either of them destroys immediately the circular trajec- 
tory. For fixed v and L, when the intensity b is not large 
enough, damping dominates and the vortex escapes from 
the system following a spiral, as explained in the previous 
Section. If b is too large, the vortex will also escape due to 
an effective drift force caused by the field, which changes 
its direction slowly enough, relative to the movement of 
the vortex. This is the case when the frequency is very 
small, such that the field is practically static. If both the 
intensity and the frequency are too large, the field will 
destroy the excitation creating many spin waves and also 
new vortices can be generated from the boundary. Many 
seemingly chaotic trajectories can be observed for high 
values of field parameters. To determine the limit cycle, 
the value of the damping is not as critical as the field 
parameters. For example, increasing the damping up to 
five times its value (e = 0.002 to 0.01) did not change 
significantly the limit cycle shown in Fig. |3 but only ac- 
celerated the reaching of it. At this point note that there 
is no resonant absorption of the energy in the ac field 
unlike the predictions of Ref. |U The field just drives 
the vortex with the frequency fi, which is always lower 
than the frequency u> of the ac field. 

All these extreme cases constrain the size and shape 
of the regimes where circular limit-trajectories appear in 
the space of field parameters (v, b). In Fig. 0] we show for 
a system radius L = 36a this parameter plane and point 
out where the vortex escapes or gives rise to a limit cycle 
or confined orbit. Similarly to what we found in the study 
of switching^ we also find "windows" , i. e. events which 
are not expected in a particular region (for instance, the 
point [y — 0.1,6 = 0.02) in the diagram). The zoom-in 
of any region of the diagram containing windows shows 
again a similar behavior. We can also observe that the 
vortex is sensitive to small variations of the field parame- 
ters, and that its behavior is not monotonous (follow for 
example the line b = 0.025 for increasing frequencies). 

When L is varied, there can appear "windows" where 
there is no limit cycle. For example for L = 36a, v = 0.1, 
b = 0.02 the vortex escapes from the system, while for 
L = 24a, 30a, on one side and L = 42a, 48a, 54a, • • • , on 
the other side, the vortex reaches a limit cycle. 

In the rest of this work we will concentrate on the 
circular limit cycle. Figs. and show the dependence 
of the vortex radial position R and azimuthal frequency 
CI as a function of the system radius L for a fixed field 
frequency v = 0.094 and four amplitudes 6. The linear 
dependence of R on the system size L is very clear from 
Fig- El for the whole range 11 < L < 56. Fig. shows 
the frequency CI of the vortex orbit as a function of l/L. 
The dependence is linear for L > 30 but not for smaller 
L indicating a possible size effect. The points missing in 
the two figures for L — 20 and b = 0.0187 correspond to 
a vortex escaping from the system. 

For a fixed system size L the features of the limit cycle 
depend on the values of the field frequency v and am- 
plitude 6. In Fig. [7] we plot the radius R of the limit 
cycle as a function of the inverse 1/v of the frequency of 
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FIG. 5: (Color online) Radius of the vortex orbit R vs. the 
system radius L, in the circular limit cycle, for a fixed field 
frequency v = 0.094 and several amplitudes b. The lines are 
there to guide the eye. Here and in the next figures R and L 
are given in units of Iq. 
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FIG. 6: (Color online) Frequency of the vortex orbit vs. 
the inverse system radius 1/L, in the circular limit cycle, for 
a fixed field frequency v = 0.094 and several amplitudes b. 

the applied field. For large frequencies one can see that 
the radius tends to a constant which is proportional to 
the amplitude b. For low frequencies the radius increases 
sharply. In this case damping plays a larger role than 
mentioned above. 

In Fig.|S|we plot the frequency f2 of the orbital motion 
of the vortex as a function of the field frequency v for four 
values of the field amplitude b. The diagonal is shown on 
the upper left corner of the picture and indicates that 

n < v. 

Although most trajectories which converge to limit cy- 
cles end up in a circular orbit around the center of the 
system, we have observed a few cases of a limit cycle 
that is not circular as shown in Fig. a) . Some chaotic 
confined trajectories can also be found as shown in the 
bottom panel of Fig. b). 
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FIG. 7: (Color online) Radius of the vortex orbit R vs. the 
inverse 1/v of the frequency v of the rotating magnetic field 
for four different amplitudes of the field. The radius of the 
system is L — 36a ~ 20. 
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FIG. 8: (Color online) Frequency of the vortex orbit fl vs. 
the frequency v of the rotating magnetic field for four different 
amplitudes of the field. The diagonal is plotted in the upper 
left corner. The radius of the system is L — 36a m 20. 



V. THEORETICAL DESCRIPTION OF THE 
VORTEX MOTION WITH ROTATING FIELD 

To describe analytically the observed vortex dynamics, 
a standard procedure is to derive Thiele-like equations, 
as it was done in Sec. IIIII without field. Due to the field 
there appears the following Zeeman term in the total 
energy (see Appendix IUl for details): 

V(r) « nbRL cos ($ - vt) . (30) 

When the vortex reaches the limit cycle, the total energy 
is constant. We have checked this fact in our simulations, 
calculating the power-dissipation relation l|A2(l . For a 
vortex, which moves according to the Thiele Ansatz, the 
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FIG. 9: Two different kinds of confined vortex trajectories 
that are not circular, occurring for large field amplitudes and 
frequencies. In the top panel the radial position R(t) of the 
vortex is periodic while it is chaotic in the bottom panel 



power-dissipation relation 1)16(1 takes the form: 



-TtrjR 2 + nbvRL sin (<I> - vt) 



The energy can tend to a constant value only when 
$ = is, so the frequency of the vortex motion should be 
equal to the driving frequency. Thus the standard Thiclc 
approach cannot provide the circular motion of the vortex 
with the orbit frequency f2 < v we have observed in our 
simulations, see previous section. The reason is that the 
field excites low-frequency quasi-Goldstone modes 30,33 , 
which can couple with the translation mode 3 ^. Therefore 
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FIG. 10: (Color online) Two trajectories of a vortex from the 
collective variable Eqs. I|37|l . starting from different initial po- 
sitions. Red line: R(0) = a, $(0) = 5tt/4, *(0) = tt/4. Blue 
line: R(0) = 10a, $(0) = 5tt/4, *(0) = it/ 4. Other parame- 
ters: b = 0.002, v = 0.125, e = 0.01, and 5 = 0.08. System 
radius: L = 78a w 44. 



it is not correct to describe the vortex as a rigid parti- 
cle and it is necessary to take into account the internal 
vortex structure. 

To describe the approach to the limit cycle we now gen- 
eralize the collective variable theory to take into account 
an internal degree of freedom of the vortex. Because the 
magnetic field changes the z component of the magneti- 
zation and generates a new ground state, it is natural to 
include into the m field an additional degree of freedom. 
To comply with the new ground state (|15f) we add to the 
</>-field (|20b|) a time dependent phase ^(r) describing ho- 
mogeneous spin precession, ^(r) can be understood as 
the generalization of an arbitrary constant phase which 
could be added in Eq. I|20b|) without changing the dy- 
namics. However, this constant phase does influence the 
dynamics if there is a constant in-plane magnetic field, 
which breaks the rotational symmetry in the xy-plane^. 

The New Ansatz that we choose is 



m(z, t) = cos 9 



\z-Z(t)\ 
I(r) 



(31a) 



[Z,T 



arg (z - Z(r)j - arg (z - Zj(t)^ 

+ argZ(r) + *(r), (31b) 

which describes a mobile vortex structure like J20J, but 
including a precession of the spins as a whole, through a 
time-dependent phase ^(r) and a dynamics of the vortex 
core, through the core width 1(t). The latter allows a 
variation of the z-component of the magnetization. Wc 
will see that in the Lagrangian the two variables / and 'J 
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are conjugate to each other so that one needs to introduce 
them together. 

We find it convenient to use in the following, instead 
of 1(t), the z-component of the total spin, 

M(r) = ~ Jd 2 xm(z : T) = M Z 2 (r), (32) 

which is related to the total number of "spin deviations" 
or "magnons" , bound in the vortex^ Here Mq 



POO 

Ala = 2 / pdp cos 9{p) « 2.75 
Jo 



(33) 



is related to the characteristic number of magnons bound 



in the static vortex. Note that without dissipation and 
for zero field, M is conserved. The field excites an inter- 
nal dynamics, changing the number of bound magnons 
and the total spin M. 

To construct effective equations we use the same vari- 
ational technique as in Section ITTT1 Besides the "vortex 
coordinates" {R, $}, we consider two "internal variables" 
{M, \ff} so that our set of collective variables is 



Xi = {R(t),^(t),M(t)^(t)} . 



(34) 



One can derive the effective Lagrangian of the system by 
inserting ansatz l|31() into the full Lagrangian Ijlljl. and 
calculating the integrals, see Appendix IU1 for details: 



J 



& „,+ r,2i. , L 2 -R 2 1 / M 

— = - i? 2 $ - In h - In — 

7T L 2 V M 



M 
Mn 



- bLRf (£) cos($ + * - vt). 



(35) 



In the same way one can derive an effective dissipative function 



J* 

7T 



(ie + i? 2 $ 2 ) f Cx + ~ ln(L 2 - R 2 ) - i In ^ + * 2 f L 2 - ^ 



Ml 



Ml 



(36) 



where the constants Ci and C2 are introduced in (|rjl3|) and (|C10fl , respectively. From the Euler-Lagrange equations 
(|26|l for the set of variables (|34|l we obtain finally 



R 

$ : 
M : 



ei? 
~2 



1 



. d + iln(L 2 -i? 2 )-iln — 
L 2 - i? 2 2i? V 2 V y 2 M 



i? 2 $ + * ( L 2 



2 V^o 



1 

M 



eC- 



M 
Mo 
M 



bLRf sin (Q + V-vt), 



'MMn 



(37a) 
(37b) 
(37c) 
(37d) 



To integrate numerically the differential algebraic sys- 
tem (ETTjl . one needs to solve at each step a linear system; 
we used the MAPLE software^ which includes such a 
facility. The set of Eqs. (|37|) describes the main features 
of the observed vortex dynamics, and yields the circu- 
lar limit cycle for the trajectory of the vortex center, see 
Fig. [TDJ Let us note that Eqs. (|37a|l and (|37bj) reduce 
to the Thiele equations for the coordinates (R, <J>) of the 
vortex center when M and \& are omitted and in this case 
no stable closed orbit is possible. Only including the in- 
ternal degrees of freedom (M, ^) one can obtain a stable 
limit cycle. 

In the parameter plane (v, b) shown in Fig. llll we indi- 
cate the two main types of trajectories found by numer- 
ical integrating Eqs. l(57)l. Vortex trajectories converge 
to a limit cycle only for b < v/2 (red domain). When 



the amplitude of the rotating field b lies above the crit- 
ical curve, the vortex escapes from the system along a 
spiral trajectory (blue domain). The model has no lower 
boundary for the limit cycle. However when the ampli- 
tude of the field lies below the critical curve (dashed line 
in Fig. II ip . the radius of the vortex orbit can become less 
than the lattice constant (green domain) . In this case dis- 
creteness effects are important for the spin system, so the 
model can no longer be adequate. 

In Fig. E| we show the radius of the vortex orbit R on 
the circular limit cycle as a function of the system size L, 
obtained from the numerical solution of Eqs. (|37|l . Notice 
the linear dependence R cx L similar to the one observed 
in the numerical simulations (see Fig. [5J . 

To analyze the main features of the model we simplify 
it, assuming that the vortex orbit is never close to the 
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FIG. 11: (Color online) The two types of trajectories ob- 
served in the (v, b) field parameter plane for the collective 
variable Eqs. JHZJ- The parameters are e = 0.01, 8 — 0.08 
and system radius L = 36a ~ 20. 
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FIG. 12: (Color online) Radius of the vortex orbit R vs. 
the system size L, in the circular limit cycle for the collective 
variable Eqs. I|37^ for a field frequency v — 0.06. The other 
parameters are the same as in Fig. 1111 



bL R 
R = rjR§ - — sinA + e— 



1 

1? 



R bL 
] R^2R 



cos A, 



M N = -eL 2 i> + bRL sin A - ei? 2 $, 
2M i Z> = N, 



(39a) 

(39b) 

(39c) 
(39d) 



The set of Eqs. (|39|l describes two damped periodically 
forced oscillators, described by two couples of variables, 
(R, $) and (N, Under the action of forcing these os- 
cillators can phase-lock and induce the limit cycle. The 
numerical study of Eqs. i|39|) reveals three different types 
of behaviors as a function of the field amplitude b for 
a fixed frequency v. We choose v — 0.06. For a small 
6 = 3- 10~ 4 , the phase A increases linearly with time, 
N oscillates and R increases very slowly without stabi- 
lization. When the amplitude is large such as b = 0.12, 
A tends to — tt, N becomes negative and then goes back 
to about 0, R increases indefinitely. For b = 3 • 10" 4 , N 
tends to a positive constant, A tends to tt so that terms 
in 7? balance and we have the limit cycle. One can see 
that the dynamics of the couple (N, ty) is fast with a typ- 
ical relaxation time of about 1/eL 2 while the dynamics 
of the couple (R, $) is slow and depends on the initial 
position Rq. The limit cycle is obtained for Rq < 0.6L, 
outside that range R increases indefinitely. 

When the solution of the system of Eqs. (|39|l converges 
to a limit cycle, we have 

R = N = 0, $ = ft = const, * = v - O. (40) 

In that case we obtain the following three algebraic equa- 
tions 



2R(v + AQ) = feLsinA, 
eL(v - fi) = W?sinA, 
-2RQ = bL cos A, 



(41a) 
(41b) 
(41c) 



where A — C\ — 1 + In L. Extracting the sin A term from 
the first and second equation, we obtain the frequency of 
the vortex motion 



system border {R <C L) and that the total z-component 
of the spin varies weakly so that N = (M — Mq)/Mo <C 1. 
Then one can simplify the expressions for the Lagrangian 
and dissipative function where the common factor 7r has 
been omitted: 



R 2 N 2 

% = M NV - i? 2 $ + — - — - bRL cos A, (38a) 



rj[R 2 + R 2 ® 2 



L 2 4 
L 2 



-4^ -heir**, 



(38b) 



where A = <f> + \& — vr, and rj was defined in (|22p. 
The equations of motion which result from l|38|) have 



CI 



1 + AR 2 /L 2 



(42) 



We now eliminate the sine and cosine terms from 141 all 
and (|41c|) . resulting in R w bL/2Cl. Combining with 142|) 
one has 



+ ^Jv 2 - Ab 2 



(43) 



This value is smaller than the driving frequency v in ac- 
cordance with our simulations. However, it is not propor- 
tional to l/L as in the spin simulations. For the radius 
of the limit cycle we have finally 



R 



bL 



+ \Jv 2 - Ab 2 



bL 
2z/ 



(44) 
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The radius of the vortex orbit R depends linearly on the 
system size in good agreement with the results of the sim- 
ulation, see Section HVl It also bears the proportionality 
to 1/v observed in the spin dynamics. 

The range of parameters, which admits limit cycle tra- 
jectories, can be estimated from the natural condition 
R < L, which gives b < 2v. However, there exist stronger 
restrictions for the limit cycle. The solution (|43[1 is real 
(not complex) only when v 2 — Ab 2 > 0. Another limit 
for the parameters is obtained from the natural condition 
RIq > a (discreteness effects are important there). Thus 
the range of parameters, which admits the limit cycle 
trajectories can be estimated as follows: 



2a b 
— < - < x 
IqL v 



(45) 



with k = 1/VJ = 1/VCi - 1 + lnL. 

For the parameters considered in Fig. ^2>r ps 0.48 so 
that the estimate (|45|l agrees with the boundary b ps v/2 
shown in the figure. 

From the above expressions one can estimate W on the 
limit cycle as 



- Vv 1 



Ab 2 



which shows that the change in magnetization ./V = 
2Mq s & due to the internal variables is small. It is never- 
theless crucial for obtaining the limit cycle. 



VI. DISCUSSION 

Another way to understand the vortex dynamics is to 
analyze the movement of individual spins. In a set of sim- 
ulations, we recorded the components of some individual 
spins to observe their time evolution. We consider a large 
enough time so that the vortex reaches the limit cycle. 
For the Fourier spectrum of the z-component of individ- 
ual spins we have observed some peaks, which appear 
naturally with the frequency of the limit cycle f2. Every 
time the vortex passes close to the spins, the spins feel 
a lick upwards. The behavior of 4>{t) for several spins is 
shown in Fig. When the vortex has reached its limit 
cycle ie for t > 500 the spins behave differently whether 
they are inside or outside the vortex orbit. Inside, <p is 
quite regular and increases linearly with time at a rate 
given by v, with <f> (po + vr. This is shown by the three 
upper curves in Fig. 1131 for r > 500 which is the time 
taken by the vortex to settle on its orbit. Outside the 
orbit and for r > 500, the increase of 4> is more irregular 
as shown by the three lower curves in Fig. 1131 There the 
Fourier spectrum of </>(t) has a main frequency ui — 17 
together with additional peaks at u ± nQ where n is an 
integer. 

Our collective variable theory describes this effect as 
we show now. We assume that the vortex has reached 
the limit cycle so that the variables <& and \& fulfil the 
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FIG. 13: (Color online) Time evolution of the c/i-field for 
spins inside and outside the vortex orbit, once the vortex has 
reached a circular limit cycle. The parameters are v — 0.1, 
b = 0.02, and L = 48a ~ 27. Inside spins ares located at 
(25a, 25a), (27a, 25a), (29a, 25a), and the outside spins are at 
(31a, 25a), (33a, 25a), (47a, 25a). 



relations 140fl . According to the Ansatz (|31b|) . on the 
limit cycle the dynamical variable <fi can be written as 

c/)(z, t) = y +VT+axg (z-Z{t)) -arg (z-Zj(t)) . (46) 

We consider the vortex to be far from the boundary, i.e. 
]J<L. Then the radius of the image-vortex trajectory 
is Rj = L 2 /R L, and for any |z| < L the last term in 
Eq. lO arg(z - Z t (t)) w tt + $(r), so 



(z, t)=<po + (v- Q)t + arg (z - Z(r)j . 



(47) 



If we consider a spin, situated at a distance \z\ > R, the 
last term in Eq. (|47J1 describes only small oscillations on 
the background of the main dependence 4>(z,t) ~ ipo + 
(y — £1)t. At the same time for a spin located at \z\ < R, 
this term is decisive. Let us consider the limiting case of a 
spin situated near the center of the system. Then arg(z — 
Z(t)) k, 7r + <I>(t), and Eq. I|47|l can be simply written 
as 4>{z, t) = (po + vt. Thus, the two regimes for the 
in-plane components of the spins are well-pronounced, 
which is confirmed by our simulations, see Fig. 1131 

In a wide range of parameters the vortex moves along 
a limit circular trajectory. When the intensity of the 
ac field exceeds a critical value, b > ;w, the vortex es- 
capes through the boundary and annihilates. This pro- 
cess is important for practical applications, because vor- 
tices are known to cause hysteresis loop in magnetic 
nanostructuresii Usually static fields are considered 
in the experiments and these cause a hysteresis of the 
M X (H X ) loop, see e.g. Refs. I6I10PT7I The satura- 
tion field in the static regime to obtain a hysteresis is 
about W0/7 (in dimensionless units 6^1). In this article 
we consider an ac driving of the vortex, which causes a 
dynamical hysteresis, M x as a function of the intensity 
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of the ac field b. Typical fields for vortex annihilation, 
b ~ nv <C 1, are much weaker than in the static regime. 
It is then much easier to destabilize the vortex with an 
ac field than with a dc field. 

Let us make some estimates. We choose permalloy (Py, 
Ni$oFe2o) magnetic nanodotaSiii. The measured value 
of M s = ~fSL 2 /a 2 = 770 G, the exchange constant A = 
JS 2 = 1.3 x 1CT 6 erg/cm, and 7 /27r = 2.95GHz/kOe. 9 
Typical fields of the vortex annihilation b ~ xv, which is 
about some tens of Oe. 

Another important fact can be seen from Fig. [21 the 
vortex is unstable in small magnetic dots, the typical 
minimal size L mnl ~ 5. For the Py magnetic dot with the 
magnetic length lg = 5.9 nm^ the minimal size for the 
vortex state magnetic dot under weak ac driving is about 
Lminlo ~ 30 ran. This means that for magnetic dots with 
diameters less than 60 nm the vortex state is unstable 
against the ac field giving rise to a single-domain state. 

In conclusion, we developed a new collective variables 
approach which describes the vortex dynamics under a 
periodic driving, taking into account internal degrees of 
freedom. To our knowledge, it is the first time that an 
interplay between internal and external degrees of free- 
dom, giving raise to the existence of stable trajectories, 
is observed in the case of 2D magnetic structures. This 
ansatz gives (up to a factor of 2) the radius of the limit 
cycle. Also the dependencies of R on the system size L, 
the field amplitude, and the frequency are correct. How- 
ever, the dependence of the vortex orbit frequency f2 on 
the system size is different from the one in the spin dy- 
namics. Moreover, in the collective variable theory the 
magnetization and vortex position variables vary on very 
different time scales, this is not the case for the spin 
dynamics. Despite this we think that this collective vari- 
able approach is very general and can be employed for 
the self-consistent description of the dynamics of differ- 
ent 2D nonlinear excitations, e.g. topological solitons in 
2D easy-axis magnets^. 



APPENDIX A: DISCRETE SPIN DYNAMICS 

While Eqs. Q are convenient for analytical con- 
sideration the presence of the time derivative on both 
sides makes them inconvenient for numerical simulations. 
Equivalent equations are obtained by forming the cross 
product of © with S n and subtracting the result from 
(J3J . In this way we get 

rf/9 f 
(1 + e 2 )^ = [S n x F n ] - | [S n x [S n x F n }} , (Al) 

where F„ — —dJ%? /dS n is the total effective field; the 
factor (1 + e 2 ) is usually neglected, or absorbed into , 
giving effective constants J, K and B. 

From the discrete dynamics IjAlf) one easily derives 
the power-dissipation relation for the total energy 34? — 
-E„S„-f„.Wehave 



dje 



E^ 



dB 
~~dt 



E F - 



dS n 
~dt 



dB 

~dt 



{l + e 2 )S 



Fn ■ [S n X [S n X 



and finally 

dJ? _ 
~~dT 



dB 



{l + e 2 )S 



E^x^-E^'^T/- (A2) 



While the first term is always negative, it is the second 
term which can give rise to transients in the relaxation 
to equilibrium, or even the resonances, depending on the 
parameters of the time-dependent magnetic field. 



APPENDIX B: COLLECTIVE VARIABLE 
EQUATIONS WITHOUT FIELD 

It is convenient to make calculations in the reference 
frame centered on the vortex whose axes are parallel to 
the standard frame 
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x - X(t) = pcosx, y- Y(t) = psinx ■ (Bl) 

Viewed from this point, the distance to the circular bor- 
der of the system changes as a function of the azimuthal 
angle x, see Fig. ^] Every integral over the domain 
\z\ < L can then be calculated as 

J f(p, a)d 2 d = (F(a)), F{a) = 2tt J f{p, a)pdp, 

\z\<L 

where a = x — $ is given by the cosine theorem 

a(a) = —Rcos a + — R 2 sin 2 a, (B2) 
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FIG. 14: (Color online) Arrangement of angles in the mobile 
frame centered in the vortex with the coordinates Z — X + 
iY = _Rexp(i<!>). 



and the averaging means (F(a)) = ^ F(a)da. We 
also have the relations 



e x cosx + e y sinx, 
- e x sinx + e y cosx, 



en 

e<s> = — e x sin <f> + e y cos $. 



e x cos $ + e y sin $, 



(B3a) 
(B3b) 
(B3c) 
(B3d) 



Then the gyroterm in the Lagrangian Sf gives 

<S X = - J d 2 ^ = -27rR.(l<j(a) - L/3] sbxa 

+ 2nR<i>([a{a) - L/3] cos a 
% = / d 2 ^m(j) — /co-R(sina) — fco -Mucosa), 



where the constant ko — 2ir J °° cos 9(p)dp. After averag- 
ing with account of the expressions 



a(a) sin a^j — 0, cos 



R 



we obtain the gyroterm in the form 



0. 



(B6) 



(B7) 



and finally, <S = -7ri? 2 $. 

Let us calculate an effective dissipative function, start- 
ing from the "microscopic" dissipative function 112fl . 
which we cut into two terms, & = &\ + ^2 with 



* = 5 fa T^- 

2 1 — m z 



h = - I d 2 i (1 - m 2 



The time derivative of the m-field can be easily calcu- 
lated in the moving frame 1|B1|) , using the traveling wave 
Ansatz (gUJl 



In order to derive an effective Lagrangian we start with 
the "microscopic" Lagrangian l|ll|) . 



int. 



d 2 i{\ - m) 



(B4) 



We will provide all the calculations for the vortex with 
unit vorticity q = 1 and positive polarity p = 1. Using 
the traveling wave Ansatz Ij20b() in the form, 

one can calculate the time derivatives in the moving 
frame (jBlf) : 



. R . R<S> 
X = — sm a cos a. 

P P 



Xi 



L 2 pR . L 2 <5> 
■ sm a 1 



R 2 pj 



Rpj 



j v p 2 i~ p 2 sin2 a - 



Here pi = \z — Zj(t)\, xi = ar g i z ~ %i ( T ) ) • In the main 
approach for R/L <C 1, one can simplify an expression 
for xii so finally we have 

4> = [k sin a - i?$ cos a) Q + . (B5) 



rn 



9' sin (^R cos a + i?<E> sin a 



(B8) 



Calculating integrals for &\ with account of (|B8f) . we 
derive: 



£7T . 



(B9) 



where A;i = 9' 2 {p)pdp. In the same way we can derive 
^2, taking into account cb from Eq. (|B5|) . 



J? 2 = e7 r |^i? 2 ^ (fc 2 + In a) sin 2 a 
+ i? 2 $ 2 ^(/c 2 +lncr) cos 2 a 
- RM ( (fc 2 + In a) sin 2a 



1 sin 2 
Using the averages 



dp - 



00 cos 2 0(p) 



sin a \n<j{a)j = ^cos a lncr(a) 
- 1 -(\na(a)) = 1 l \n(L 2 -R 2 ) 



dp. 



(BIO) 



(Bll) 
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we calculate the dissipative function in the form 



£7T 



C\ + ^HL 2 -R 2 ) 



(tf+R 2 ® 2 ). (B12) 



Here the constant C\ = ki + ki, 

„ 5 f 1 sin 2 6{p) , 

Ci = - + / K -^dp - 

4 Jo P M 

6' 2 {p)pdptt 2.31. 



cos 2 



dp 



(B13) 



Supposing that the vortex is not close to the boundary, 
i.e. R <C L, we obtain the effective dissipative function 
in the form (|25p. 



Besides this direct influence on the system, the mag- 
netic field also changes the gyroterm in the effective La- 
grangian, and the energy of the system. These changes 
result from the internal motion of the vortex through 
£(t), and from the uniform spin precession through ^(r). 
This does not change the gyroterm Sfi, which has the 
same form as in (IB7|l , but there appears the contribution 
Sf 2 = M\fr. This can be easily calculated with account of 
the time derivative 

4> = i>+ (Asina- RQcosa) (j^ + J^J ■ ( C6 ) 

The total energy functional © can be written in the form 
S = Sx + S 2 + S-i + V with 



APPENDIX C: COLLECTIVE VARIABLE 
EQUATIONS WITH FIELD 

First we calculate an effective Zeeman energy for the 
standard Thiele-like motion of the vortex. Inserting the 
traveling wave Ansatz Il2()[l into the "microscopic" Zee- 
man energy @), and calculating the integrals, we get the 
effective energy in the form: 



V(r) 



1 f 27T 

2 b J d x[^ 2 (x-*)-ci]cos( 



wbRLf(2) cos ($-z/r) 



where 



(CI) 



(C2) 



6i = - a £- - = km, 



1 — m 2 



(C7a) 



\ J ^ (1 _ m 2 )(v ^ )2Rj7rln ^_^! ; (C71)) 



i(r)L 



1 r -2, 2 



(C7c) 



The term ^2, which describes the interaction between the 
vortex and its image, can be derived from (|23|l . simply re- 
placing Zo by Z(t). In the last anisotropy term ^3 we have 
used the relation J" 00 cos 8 2 (p)pdp = 1/2, see Ref. E2I 
Combining all terms of the Lagrangian and omitting the 
constant term <§\ , one obtains the effective Lagrangian of 
the system (1531 . 

The dissipative function contains two dynamical con- 
tributions. The first one is due to the time dependence 
of the m-field: 



where E(x) and K{x) are elliptical integrals. When the 
vortex is far from the boundary, which is the case of 
interest, one can expand this function into the series, 



1 



In the main approach it leads to the 



Zeeman term in the form H30|) . The corresponding mag- 
netic force 



F h = -V R V = e R nbLg (f ) cos ($ - vr) 



e^wbRLf (£) sin ($ - i/r) 



(C3) 



where the function 

g(x) = f{x) + xf'(x) 



4 

3^ 



K(x) 



3-1) 



(C4) 



1 - 



For x <C 1 it has the following expansion g(x) 
3x 2 /8. 

Let us calculate the same Zeeman energy using the new 
Ansatz <|31|) . One can derive a Zeeman term similar to 

EH 



V{t) = nbRLf (f ) cos ($ + * - i/r) 



(C5) 



0' sin 6» / M 



— — p + Rcosa + R$sma . (C8) 



f(r) V2M' 

This term ^ can be derived in way similar to 1B9|I : 



1 f°° 

C 2 = - / e' 2 (p)p 3 dp sa 0.48. 

2 Jo 



(C9) 
(CIO) 



To calculate the second term J^2 we use from Eq. (|C6jl 
and obtain 



57T 



jr 2 = r-<; (i? z + i? 2 $ 2 



1, L 2 -R 2 

k2 + 2 ln n^r 



(Cll) 



^ 2 [L 2 -Z 2 (r)] +2 J R 2 $<j 



The total effective dissipative function & = &\ + has 
the form |J5SJ|. 
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